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We theoretically analyse the possibility to electrostati¬ 
cally confine electrons in circular quantum dot arrays, 
impressed on contacted graphene nanoribbons by top 
gates. Utilising exact numerical techniques, we com¬ 
pute the scattering efficiency of a single dot and demon¬ 
strate that for small-sized scatterers the cross-sections 
are dominated by quantum effects, where resonant scat¬ 
tering leads to a series of quasi-bound dot states. Calcu¬ 
lating the conductance and the local density of states for 
quantum dot superlattices we show that the resonant car¬ 
rier transport through such graphene-based nanostruc¬ 
tures can be easily tuned by varying the gate voltage. 



Schematic representation of a Dirac electron wave packet im¬ 
pinging on a circular, electrostatically defined quantum dot. 
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1 Introduction High-quality graphene nanostruc¬ 
tures are the most likely building blocks in future electron¬ 
ics, plasmonics and photonics assemblies. Most of their 
striking features arise from the strictly two-dimensional, 
honeycomb lattice structure of the basic material, which 
causes the nontrivial topology of the electronic wave func¬ 
tion and an almost linear low-energy spectrum of the 
chiral quasiparticles (charge carriers) near the so-called 
Dirac nodal points m. From an application-technological 
point of view, the tunability of the transport properties 
of graphene by external electric and magnetic fields is 
of particular importance O. This allows a controlled 
modification of selected areas of the sample by gating. 
For instance, applying nanoscale top gates on graphene 
nanoribbons (GNRs) or bilayer graphene, single or double 
quantum dots with particular shape have been produced in 
recent experiments ma. 


For a large circular dots, the refraction at the boundary 
gives rise to two coalescing caustics that focus the electron 
density in the disk-shaped region I6ll7]| . Resonances in the 
conductance m and the scattering cross section 13 indi¬ 
cate electron confinement for small quantum dots in mono- 
layer graphene as well. Apparently the difficulty to localise 
Dirac electrons by electrostatic potential barriers (which 
is related to the Klein tunnelling phenomenon) could be 
overcome: Because of the circular symmetry, the angle 
of incidence at the dot’s boundary is a constant of mo¬ 
tion so that electrons with nonzero angular momentum are 
confined P8(. Forward scattering and Klein tunnelling will 
also be suppressed if a Fano resonance arises between the 
background partition and the resonant contribution to elec¬ 
tron scattering m. These results, obtained within contin¬ 
uum Dirac theory, were numerically confirmed for a lat¬ 
tice model mm . The role of resonances have been anal¬ 
ysed in various graphene scattering experiments II121I13L 
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Note that the scattering (transport and optics) of massless 
Dirac fermions on potential barriers substantially differs 
from that of massive chiral Dirac fermions or Schrodinger 
electrons m. 

2 Model and numerical methods In this work, we 
address the scattering and (linear) transport in armchair 
GNRs with gate-defined quantum dots in the end-contacted 
lead-sample geometry most relevant for experiments (TSl 
fTMTTt . The leads can be viewed as semi-infinite quantum 
wires; the contacts are metallic or gated graphene. 

We describe the sample by the tight-binding Hamilto¬ 
nian 

F = y]lAclc,-i^(clc^.+H.c.), (1) 

* (ij) 

where annihilates (creates) an electron at Wannier 
site i of a honeycomb lattice with carbon-carbon distance 
a = 1.42 A. Here the hopping amplitude between near¬ 
est neighbour sites {ij) is given by t (:^ 3eV) and the 
site-dependent potentials Vi = V^^0{R — |r^ — r^|) 
model electrostatically defined quantum dots with radius 
R centred at positions r^. 

In the limit of vanishing bias voltage, the conductance 
between the left (L) and right (R) leads can be obtained 
within the limits of the Landauer-Biittiker approach: 

GlR = Y ^ |5'n,mP , (2) 

nEL,mER 

where Sn,m is the scattering matrix between all open (i.e., 
active) lead channels m . We furthermore exploit the local 
density of states (LDOS) at a certain site i of the sample, 

LDOS(E) = (3) 

I 

as a probe of a possible electrostatic confinement of the 
charge carrier I191I2QI . In Eq. ([^, E is the particle’s en¬ 
ergy and the sum extends over all single-electron eigen¬ 
states |/) = c]^ |0) of H with energy Ei. We use the kernel 
polynomial method ED and the ‘Kwant’ software pack¬ 
age ^T2\ for evaluating the conductance and LDOS numer¬ 
ically. 

3 Results and discussion 

3.1 Single quantum dot We first consider a single 
dot, see Fig. To begin with we determine the scattering 
efficiency Q, which is the scattering cross section divided 
by the geometric cross section, of a plane Dirac electron 
wave hitting on a dot embedded in an infinite graphene 
sheet, within the Dirac approximation (for details of such 
a calculation cf. Ref. 13). Since the energy is conserved 
in our setup (i.e., there is no coupling to bath degrees of 
freedom etc.), the scattering is always elastic (coherent). 
The data displayed in the upper panel for = 0.051 and 
R = b nm shows that a series of resonances appears when 
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Figure 1 Total and partial scattering efficiencies (top), 
conductance (middle) and LDOS (bottom) of a single 
quantum dot (the black circle specifies the dimension of 
the dot). For further explanation see text. 


the gate voltage is varied (besides the rather broad uq and 
ai modes, the a^>i resonances are very sharp), which 
can be attributed to the excitations of higher dot normal 
modes IQ. So each maximum marks a combination of E, 
R, and V, where a certain normal mode fits particularly 
well into the dot, while all other modes are suppressed. The 
middle panel indicates that these signatures also show up in 
the conductance, which this time, however, was computed 
for a 20 nm wide armchair GNR placed between leads. 
That is why we can assign the signal at ^ — O.OQt to 

the ai mode. The value for G[V = 0] tells us that the sys¬ 
tem exhibits three open transport channels. The lowermost 
panel presents the LDOS close to the ai resonance marked 
in the G ys. V plot by the blue diamond. Obviously the 
electron is almost completely localised at the quantum dot. 
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Figure 2 Conductance (top) and LDOS (lower panels) of a 
GNR with a linear array of quantum dots. The colour code 
is the same as in Fig.[^ 


where the superposition of the states with j = ±(m +1/2) 
leads to the characteristic vortex pattern of the ai mode. 
There are six vortices close to the boundary of the dot 
which are reflected by three preferred scattering directions 
in the far field (91. 

3.2 Qunatum dot chain We now study quantum dots 
arranged like a string of pearls. Figure [^presents the con¬ 
ductance and the LDOS for a contacted GNR sample with 
five dots at intervals of 20 nm. Again we set = 5 nm, 
E = 0.05t, and the width of the armchair GNR is 20 nm. 
Note that the conductance undergoes a dramatic change 
if V is slightly varied near the resonance points (see red 
and blue diamonds). In this way the system may act as a 
switch. The results for the LDOS provide further insight 
into how the system behaves locally for V = —0.088 f 
(V = —0.0801), i.e., at the local minimum (maximum) of 
the conductance. At 12 = —0.0881, E = 0.05f (middle 
panel), the quantum dots are off-resonance from an eigen¬ 
mode of the particle: As a result the scattering is enhanced 
and the conductance is low (note that the figure shows a 
stationary situation). For an applied bias 12 = —0.0801 
(lower panel) all dots largely support the ai mode. Since 
the wave functions of quasi-bound ai modes at neighbour¬ 
ing dots overlap resonantly, electrons are no longer perma¬ 
nently confined to a particular dot but will be transferred 
through the GNR by an effective inter-dot hopping process. 
Remarkably, in this case the conductance is almost as large 
as for the GNR without quantum dots (G[12 = 0]) where 
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Figure 3 Conductance (top) and LDOS (middle/bottom) 
of an armchair GNR with a square dot superlattice. 

the LDOS is the same (i.e., uniform) everywhere between 
the leads. 

3.3 Quantum dot matrix Figurej^illustrates that the 
same behaviour could be achieved for a square quantum 
dot superlattice made on an armchair GNR (now the width 
of the GNR is 60 nm, the other parameters are the same 
as in Fig. [^. For 12 = 0 the system has nine open trans¬ 
port channels. Again the relatively large conductance at 
12 = —0.0801 (compared to G for 12 = —0.0651) orig¬ 
inates from an effective (coherent) inter-dot transfer of the 
electron. This hopping process is expected to take place on 
a reduced time scale CD. As can be seen from the lower¬ 
most panel, at 12 = —0.065 f, the LDOS notably shrinks if 
we move to the centre of the GNR; this depletion of the lo¬ 
cal particle density reduces the inter-dot hopping and con¬ 
sequently the conductance. 
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Figure 4 LDOS for a linear (top) and a square (bottom) 
quantum dot array on a GNR with almost vanishing con¬ 
ductance. 


While graphene has a gapless band structure and also 
the armchair GNRs considered in Figs. [T][^ have G > 0 for 
all E at V = 0 (i.e., there are always open transport chan¬ 
nels), the conductance of the GNRs with imprinted quan¬ 
tum dots may tend to zero in certain ranges of V ; see the 
black squares in Figs. and Apparently the linear and 
square lattice periodic arrays of gate-defined quantum dots 
give rise to new superlattice band structures, which exhibit 
energy (pseudo) gaps just as for ordinary solids 1^ . This 
is reflected in the LDOS depicted in Fig.[^ 

4 Conclusions To summarise, using exact numerics, 
we have studied the electronic and transport properties of 
contacted graphene nanoribbons with gate-defined quan¬ 
tum dots of realistic size (i.e., with a few thousand carbon 
atoms that feel the electrostatic potential step). The con¬ 
ductivity of such systems is very sensitive to the model 
parameters and can even vanish. This provides an oppor¬ 
tunity to use them as switches by varying the top gate po¬ 
tential. A detailed analysis of the local density of states 
shows that an electrostatic electron confinement is possi¬ 
ble for single quantum dots in the quantum regime due 
to resonant scattering. Thereby the incident electron is fed 
into vortices, which trap the particle. That means for small 
dots the particle confinement is not by total internal refiec- 
tion. In addition, for periodic quantum dot arrays a super¬ 
lattice band structure evolves that allows for an effective 
inter-dot transport on a reduced energy scale. The phe¬ 
nomena detected in this work might be exploited in ac¬ 
tual graphene-based devices. From a theoretical point of 
view, it seems worthwhile to study a modulation of the 
quantum dot barriers—via the top gates—by oscillating 


fields, which induces both inelastic tunnelling and reso¬ 
nant transport through excited states ll^ . Another issue is 
the electrical and magneto transport trough quantums dots 
on folded graphene samples that exhibit superlattice struc¬ 
tures in itself flSj . We finally like to emphasize that such 
graphene quantum dots might used to construct logic gates 
and “quantum computing” devices (261, which would be 
also an interesting topic for future studies. 
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